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Abstract 

The issue of why computational resolution in Navier-Stokes turbulence is so hard to achieve 
is addressed. It is shown that Navier-Stokes solutions can potentially behave differently in two 
distinct regions of space-time where is comprised of a union of disjoint space-time 
'anomalies'. Large values of |Va;| dominate K~, which is consistent with the formation of vortex 
sheets or tightly-coiled filaments. The local number of degrees of freedom needed to resolve 
the regions in satisfies 

where TZu ~ uL/v is a Reynolds number dependent on the local velocity field u(a;, t). 
1 Introduction 

The space-time distribution and morphology of the vorticity and strain fields in three-dimensional 
Navier-Stokes turbulence has remained a puzzle since Batchelor and Townsend discovered exper- 
imentally the phenomenon of intermittency in experimental flows. Instead of observing Gaussian 
behaviour in the flatness factor and similar quantities, they discovered the spiky spectra which are 
now recognized as typical for intermittent turbulent flows |21 151141171 [51 191 [TUl [IT] . 




Figure 1: The picture, courtesy of Jorg Schumacher of the Technical University llmenau, is a three-dimensional statistically 
stationary homogeneous isotropic flow at a Taylor microscale Reynolds number of 107 showing iso-surfaces of at a level of 10 
times the average of lS^ . The cube with side-length 27r is resolved with 2048'^ grid points which translates into 3 grid spacings 
for the Kolmogorov length A^. 

The application of colour graphics in this past generation has dramatically illustrated how the mor- 
phologies of the vorticity and strain fields are typically dominated by 'thin sets' : see [T21[T5]- These 
usually form initially as quasi-two-dimensional vortex sheets which, under interaction, roll up into 
a tangle of quasi-one-dimensional tubes [14]. It is also important to note that vorticity and strain 
accumulate on significantly different sets [T5]: indeed, there has been some debate over their relative 



importance [16l [TTl, 17]. Figure 1, courtesy of Jorg Schumacher, is a snapshot illustration of the vor- 
ticity (enstrophy) field . Experiments show that these structures spontaneously appear and disappear 
[HEKH]. While there exists an extensive literature on intermittency in the statistical physics literature 
concerning Kolmogorov's theory (^KTUlttl], no satisfactory theoretical explanation for the high de- 
gree of space-time complexity of these phenomenon has ever been given based on three-dimensional 
Navier-Stokes solutions, nor has any mathematical explanation been forthcoming why vortex sheets 
should be, at least initially, the favoured topology. The consequences of this behaviour is far-reaching. 
The spontaneous appearance of these structures, often at very short length scales, creates severe 
resolution problems which have not been been wholly solved despite the increase of computing power 
in this past generation. To provide a partial theoretical explanation of this is one of the main aims 
of this paper. Landau's heuristic estimate for the number of degrees of freedom M ~ Re^^^ needed 
to resolve a turbulent flow is based on space-time averages and is by no means enough to resolve 
the thin structures discussed above: for instance see Kerr [15j, Schumacher, Sreenivasan ii Yakhot 
[20] , Schumacher, Sreenivasan ii Yeung [21] and Sreenivasan [22| . 

Turning now to the status of Navier-Stokes solutions, there are generally two prevailing views. 
The first, which is generally held by the computational fluid dynamics community, is that the Navier- 
Stokes equations are regular: that is, it is believed that unique solutions exist that can ultimately 
be resolved provided enough computing power is made available in the future. This is equivalent to 
the assumption that the Navier-Stokes equations possess strong solutions. The second view, held 
more by Navier-Stokes analysts, is that the unsolved regularity problem leaves open the possibility 
of singularity formation [231 [211 ES [21] ■ Caffarelli, Kohn and Nirenberg [27] have shown that the 
potentially singular set has zero one-dimensional Hausdorff measure, which means that if singularities 
do occur in space-time then they must be rare events. Various routes to the efficient construction 
of suitable weak solutions are in evidence f2Sl 1221 EHl [S]- 

For the purposes of providing a mathematical explanation for the resolution problem outlined 
above, the first view will be taken in this paper. In fact it cannot be emphasized enough that a flow 
may be regular but could nevertheless be highly intermittent thus rendering the singular set empty. 
To make these arguments accessible to a wide readership the proofs, although minimal, are relegated 
to appendices which can be ignored if the reader so wishes. 



2 Results based on space-time averages 

The setting is the following: we consider the incompressible {dWu = 0), three-dimensional Navier- 
Stokes equations for the velocity field u{x, t) with mean-zero, divergence-free forcing, 

ut + u - Vu = uAu -Vp + f{x) , (2.1) 

on a periodic three-dimensional domain V = [0, L]^. Leray's energy inequality [25l [26] 



gives some information on the energy ||m||2- In (2.2) || • ||2 is the L -norm on the periodic volume 
V = [0, L]^. The method of Doering and Foias [32J is now applicable in which the space-time 
averaged velocity U and the energy dissipation rate e are found to be a priori bounded quantities 

U^ = L-m\u\\l) e = uL-mMl), (2.3) 



2 



in which the symbol ( • ) for the long-time average is 

1 /■* 

{g) =limsup- / g{T)dT. 



(2.4) 



For forcing concentrated around one length scale £ with £ taken for simplicity to be £ = L/2t:, the 
Reynolds and Grashof numbers are defined by 



Re 



U£ 



Gr 



£ /r 



in which fr,ns = L ^/2||/||2. (2.2) gives 



e < u^L-^GrRe . 



(2.5) 



(2.6) 



Moreover, in [32J it has been shown that at high values of Gr Navier-Stokes solutions obey Gr < Re^ 



and so the right hand side of (2.6) can be estimated as GrRe < cRe^. Because U and e are bounded 



quantities, the definitions of the respective inverse Taylor micro-scale and Kolmogorov lengths A^,^ 



and A, ^ are on a sound footing: 



A;; 



(ll-lli) 



u 



Then it is easily shown that [32| 

LKj,^ < cRe^/'^ , 



)' 



la; 



1/4 



(2.7) 



(2.^ 



This leads to an estimate for the number of degree of freedom in the system J\f{Ak) <cRe^/^ based 
on the number of small vortices of volume A| relative to the box volume . This is consistent with 
Landau's heuristic and Kolmogorov's scaling arguments [TUl [TT]. These space-time averages are 
extremely useful in setting average magnitudes relative to which other quantities can be measured 
but they hide strong spiky variations in local behaviour. Two new results, which are based on an 
extension of (2.6) to higher moments (m > 1), are proved in Appendix|A] 




and 




cjP'" dV 



Vc^P'" dV 



< CQ^rnV'^^-^L'^Re'^ 



(2.9) 



(2.10) 



In (2.10) when m = 1 the exponent on the integral within the time average is not unity but 1/3. 



These averages hint at some control over large fluctuations in u but not enough information is 



available to understand the behaviour of local space-time variations. Nevertheless, (2.10) plays an 
important role and will be referred to later in O 
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3 Local space-time results 



The task is now to consider how a fluid can behave in local regions of space-time based on the 
assumption, discussed in ^ that strong solutions exist. This approach requires some differences of 
definition, particularly in the Reynolds and Grashof numbers of the last section. Re and Gr, whose 
definitions were based on a space-time average in the case of Re and a spatial average in the case 
of Gr. Now we define local Reynolds and Grashof numbers as 

„ , , L\u(x,t)\ ^, , L^lffa;)! 

The idea is to consider the enstrophy and palenstrophy 



Hx{t)= / \uYdV H2{t)= / \Vu:\^dV. (3.2) 

Jv Jv 

and derive, as in Appendix[B] a differential inequality involving these (valid for 6 > g > 3) 

This is a modification of a version first proved by Ladyzhenskaya [24] in a regularity proof conditional 
upon the velocity u being bounded in \\u\\q for g > 3. Instead, when q = Q, the nonlinear term is 
proportional to the perfect integral ||it||g. Time averaging as in (2.4) and re-scaling so that each 
term is dimensionless, (3.3) converts to a 4-integral 

C67^^-£l|VcJ|2 + 2C?2| >0. (3.4) 

A trivial but important observation is that a positive integral need not necessarily have an integrand 
that is positive everywhere in space-time. Some consequences of this are: 

1. There are regions of space-time R+ C M"^ on which 

^2 

<CQnl + 2g\ (3.5) 

2. In contrast, there are (potential) disjoint regions of, or 'anomalies', in space-time on which 

-^\Vu\^>cenl + 2G\ (3.6) 

HuJq 

The union of these anomalies we call C M^. The contribution from cannot be too 



large given the positivity of the 4-integral in (3.4). 

Given that the set is non-empty, the first observation is that the very large nonlinearityTZf^ 
amplifies the response in the magnitude of |Vcjp to relatively small fluctuations in the local 
velocity field u. Using the fact that the higher moments of |Va;| are controlled as in (2.10), it 
is clear that |V(x;| cannot become too large everywhere in or this average will be violated. 
Very large behaviour in must therefore be balanced by smaller behaviour in R+. 

The product of this amplified and highly uneven response is an intermittent spectrum. This 
is consistent with the remark of Batchelor and Townsend [IJ where they suggested that large 
wave-number components are concentrated in isolated flow regions with an uneven energy 
distribution associated with the small scale components. 
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4. The second observation is that the sudden rapid increase of the gradient \ Vu\ as one moves 
from a region in R+ across into an anomaly in M~ is consistent with the formation of vortex 
sheet-like structure^ or perhaps tightly-coiled filaments. The subsequent roll-up of these 
sheets when they interact, as observed in numerical experinnents, is not explained but since 
both topologies have a small packing fraction the roll-up of one into the other would not 
change this. 

5. Defining a local Kraichnan length as 

= ^ = ^-^> (3-7) 



then (3.5) and (3.6) can be re-written as 



(LXkrf ^8cenl + lQg^ on R± . (3.8) 
Thus, to resolve the structures in R^, the number of degrees of freedomAA^ is estimated as 

Af^^c±{nl + g) on R± (3.9) 
compared with the standard Re^^^ needed on the space-time average. Moreover, TZu in par- 



ticular, might be considerably larger than Re in or near some anomalies. Given that (3.9) is a 
lower bound in R^, might be considerably larger than TZ^, showing why resolution might 
be lost locally, if u becomes too large. 

4 Summary and discussion : vorticity versus strain 

The arguments of the previous section show clearly that for strong solutions of the Navier-Stokes 
equations, space-time can potentially be split into two parts R^. While the closed set R+ is dominant 



because of the positivity of the space-time integral in (3.4), the vary lare lower bound TZ^ in (3.8) 



provokes high values of {VujI in the space-time anomalies, whose union is the open set R . The 



control of space-time averages of higher moments in (2.10) insists that large values of |Va;| in 
R~ must be balanced by smaller values in R+. The violent increase of u in the relatively small 
space-time region occupied by these anomalies suggests the formation of vortex sheet-like structures 
or tightly-coiled filaments. Numerical evidence suggests that sheets roll up into tubes when they 
interact, as in Figure 1 [Hj. Indeed, both topologies have an appropriately small packing fraction 



which is no doubt why they are commonly observed. The bounds on Af^ in (3.9) are dependent 
on the velocity field at local space-time points of the flow. Despite the assumption of regularity of 
solutions, TZu may potentially reach large values, making the 7^,^ lower bound enormous. This would 
account for local difficulties in resolution and illustrates the need to carefully monitor values of 7^^ 
in a numerical calculation : see, for instance, [20| . 

These ideas can also be used in an alternative manner to see the effect of strain : see [TBI [iTl [T^ 17] 
and [12]. In this case it is more appropriate to define different local Reynolds numbers as 

= , Tip = . (4.1) 



^Formally, a vortex sheet is considered to have formed when a jump occurs in tangential vorticity as one moves in 
a normal direction ; in our case we are unable to make any distinction between the normal and tangential directions. 
See Majda &i Bertozzi [34J for a discussion and references on the formation of vortex sheets in the Euler equations. 
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In (4.1) ps{x, t) is the spectral radius of the strain rate matrix S which appears because u-u-Vu 



uj ■ Su. The equivalent of the 4-integral in (3.4) is 



2npnl-^\vu\^ + g^\ dv)>o. (4.2) 

/V L ^0 J / 

Similar conclusions can be reached to those of f|3] regarding the effect the strain field and vorticity 
fields have on |Va;| in regions of space-time. These will be different from those regions contained 

- the 



m 



of (3.4). An estimate for the number of degrees of freedom needed to resolve 



equivalent of (3.9) - is 



(4.3) 
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A Proof of (2.9) and (2.10) in §2 



The estimates in this appendix are based on the assumption that strong solutions of the Navier- 
Stokes equations exist. Consider the result of Foias, Guillope and Temam [33] for time averages of 
higher semi-norms which is deliberately written for our present purposes in terms of Gr and Re 



Hr, 



A Sobolev inequality for m > 1 gives 



l-a 



where a = 3{m — 1)/Am. Thus, taking n = 3 in ( A.l ), we have 



2m 
I 4m-3 

I 2m 



< C 



15(m-l) m+3 
^l/5\ 4(4m-3) ^4(4m-3) 



15(m-l) 
l/5\ \ 4(4m-3) 



m + 3 



The averages of |Va;| in (2.10) can be found in the same way using n = 4 in (A.l) 



(A.l) 



(A.2) 



(A.3) 



B Proof of Ladyzhenskaya's inequality (3.3) 



Consider Hi defined in (3.2). For strong solutions 

Hi= u ■ {uAu + u; ■Vu + curl /} dV . 
Jv 

Using the divergence theorem to estimate the vortex stretching term we find 



X 



u ■ iu ■ Vu dV 



< \\Vuj\\2\\uj\\s\\u\ 



(B.l) 



(B.2) 
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where q ^ + s ^ = |. A Sobolev inequality on \\uj\\s gives 



\UJ\ 



< ||Vu||s < c||V^ti|||||ti||^ 



l-a 



(B.3) 



where a is given by a = There are restrictions on a in the form ^ < a < 1, which means that 

3 < g < 6. Thus we have 



uj ■ u ■ Vu dV 



< c||Vcj||^+°i|u||2-'^ 



, I _ 1+2: ,, ,, 1 



l-a 



l + g 



2(2-a) 



l-a 



( B.l ) becomes 



1^1 < -'^ 



3q 



-^]H2 + CgU 'i-^\\u\\^ +11/112^^2 



We also know that 



g-3 
q + 6 



.1/2 



g + 6 
g-3 



.1/2 



2(g + 6) 



2(g-3) 



Together (B.5) gives 



' 2(g + 6) ^ II 119 ^ 2(g-3) 



--^ll/lli, 



where one recalls that 6 > g > 3. This is (3.3) 
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